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Abstract 

The exact solution of the Dorokhov-Meno-Pereyra-Kumar-equation for quasi 
one-dimensional disordered conductors in the unitary symmetry class is em- 
ployed to calculate all m-point correlation functions by a generalization of the 
method of orthogonal polynomials. We obtain closed expressions for the first 
two conductance moments which are valid for the whole range of length scales 
from the metallic regime [L <C Nl) to the insulating regime {L ^ Nl) and for 
arbitrary channel number. In the limit iV — > oo (with L/{Nl) = const.) our 
expressions agree exactly with those of the non-linear a-model derived from 
microscopic Hamiltonians. 
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Typeset using REVT^ 



The current understanding of transport and localization in quasi one-dimensional dis- 
ordered wires is essentially based on two competing approaches. The first one is a 
random matrix approach in which the transfer matrix is written as a product of statistical 
independent building blocks modelling the quasi one-dimensional structure of the wire. The 
distribution of each building block is calculated by maximizing its entropy (local maximum 
entropy approach). This procedure relies on the assumption that all scattering channels are 
equivalent. The probability distribution p{{\i}, L) (the Aj denote the usual radial coordi- 
nates which are related to the transmission eigenvalues via Tj = (1 + Aj)~^) is determined 
by a Fokker-Planck equation (henceforth called DMPK-equation) describing the evolution 
with the length L of the wire : 

dX. J{{Xi}) ' ^ ^ 

>/({M) = ni^-^-A,i'' ■ (2) 

i>j 

The three possible values of /? = 1,2,4 correspond to the orthogonal, unitary or symplectic 
symmetry classes, respectively, and ^ = {PN + 2 — /5) / is the localization length. 

An alternative description of quasi one-dimensional wires starts from a more direct, mi- 
croscopic approach. The metallic regime (Z <C L -C Nl) is treated in terms of diagrammatic 
perturbation theory [§,0. To extend the results to all length scales it is necessary to con- 
sider the one-dimensional supersymmetric non-linear cr-model [p|-p!l||. Developing the Fourier 
analysis on the matrix space defining the cr-model, the average of the conductance and 
its variance have been calculated for all length scales from the metallic to the insulating 
regime (L 3> Nl). Since the one-dimensional cr-model can be derived from a Hamiltonian 



with white-noise potential P], Wegner's A^-orbital model |]lTl, and a random band matrix 



approach for the Hamiltonian |]12|, it provides a description of quasi one-dimensional disor- 



dered conductors that is rather independent of microscopic details. 

The results for the universal conductance fluctuations and weak-localization corrections 
in the metallic regime obtained by the Fokker-Planck approach and by the cr-model (and 
also by the diagrammatic perturbation theory) are in agreement. In the insulating regime, 
the behavior of the probability distribution following from the DMPK-equation is rather 
well understood. The Lyapunov exponents are known analytically and the Aj (and hence 
the conductance g = J^A^ + Aj)"^)) follow a log-normal distribution |T3HT5|]. However, 



a direct comparison with the results of the a-model is rather difficult because only the 
averages (g) and (g"^) (instead of (iTig)) are known. The issue is rather important due to 
the unusual results of the cr-model for the symplectic symmetry class (systems with time- 
reversal symmetry and rather strong spin-orbit coupling): Both the average conductance 
and its variance were found to have finite values as L/C, oo QJl^l- It is therefore of 
considerable interest to derive from the DMPK-equation closed expressions of the average 
conductance in the insulating regime, even though {g) does not coincide with the typical 
conductance exp{{\ng)). 

The situation in the unitary symmetry class (systems with no time reversal symmetry) 



is particularly favorable because recently |jT6| the DMPK-equation has been solved exactly 
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by means of a Sutherland-transformation |^ leading to an independent Fermion problem. 
Based on the expression of p{{Xi},L) given in Ref. [16|, this letter is concerned with the 
exact calculation of the m-point correlation functions 



Rm{Xi, . . . , Am, L) 

m 



(3) 



(A^ - m)\ Jo 



dX 



m+l 



■dXN p{{Xi},L) 



that are usually considered in the theory of random matrices |T8 



The exact result for the joint probability distribution given in Ref. can be expressed 
in terms of the variables A,- as 



p{{Xi},t) oc Y[{Xi - Xj) det{gra-i{Xj,t))e 
i>j 



(4) 



with 



5'm(A,t) 



rfA;(P)"^ |tanh(^) 



xPi(,,_i)(l + 2A)e 



(5) 



and a constant to be found in [|l6i|T9| . The functions Pi(j^_]^^(l + 2A) denote the gener- 
alized Legendre functions. We have introduced the abbreviation t = i.e. the system 
length measured in units of (twice) the localization length. 

In random matrix theory, one typically applies the method of orthogonal polynomials 
18| to calculate the m-point functions (^. This method cannot directly be carried over to 



the present case because the pairwise "interaction" between the A, eigenvalues is no longer 
purely logarithmic |T6|. A first generalization |^ of the method is based on the application 
of biorthogonal polynomials. Extending the idea of Ref. [2C], we write the probability 
distribution (11) product of two determinants 



p{{Xi},t) ocdet{Qn-i{Xj,t)) det{hm-iiXj,t)) 



(6) 



with "dual" functions QniX, t) and hm{X, t) (n, m = 0, 1, . . . , A^ — 1) that fulfill the biorthog- 
onality relation 



dX Qn{X,t) hmiX,t) = 6n 



(7) 



Expressing the product of the Aj-differences in (^ as a Vandermonde determinant with 
(arbitrary) polynomials Qn{X,t) of degree n, we arrive at (^). The biorthogonality (|^) can 
be achieved by a suitable linear transformation of the functions gm{X,t). In the following, 
we choose the Q„(A,t) as 



Q„(A,t) = P„(1 + 2A) e 



-ent 



(8) 



where Pn are the Legendre polynomials. 

The next step is the evaluation of the integral of each combination of Qn{X,t) with 
gm{X,t). The Legendre functions Pi(j^_]^)(l + 2A) (see Eq. (|^)) are eigenfunctions of the 
differential operator 
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with eigenvalue k"^ . The operator D is related to the one particle Hamiltonian Tio used in 
Ref. by a suitable transformation. In the following, we need three properties of gm{^, t), 
namely 

gmiKt) = D"^goiX,t) , (10) 
^^go{X,t) = -DgoiX,t) , (11) 
goiX,0) = 5iX) . (12) 

Eqs. ([I0|), { \n\ ) follow directly from (^. Eq. (|I2|) reflects the initial condition for the one 
particle Green's function ||16i. In order to proceed, we need the identity 



/ dXQn{X,t)gUX,t) = {-er,r , (13) 
Jo 

which follows from the fact that the Legendre polynomials P„(1 + 2A) are also eigenf unctions 
of the differential operator D (with eigenvalue — £«)• One can verify for m = via (|TTl , p!2[) 
that the integral in ([TB|) does not depend on t and is equal to unity. The generalization to 
arbitrary m follows directly from Eq. (0). 

We now construct the functions hm{X,t) in analogy to Eq. (^), i.e. 

hm{X,t) = r dkLm{k^) |tanh(f ) 

•J 

xPi(,,_i)(l + 2A)e-^'* . (14) 

Here, the Lm{- ■ ■) denote a set of linearly independent polynomials of (maximal) degree 
— 1. The integral (13) with the gm replaced by hm takes the value Lm{—£n)- The 
biorthogonality relation (|^) is therefore fulfilled if we choose the polynomials Lm as the 
Lagrangian interpolation polynomials 

Lm{z)= U ■ (15) 

As a key step to apply standard methods from random matrix theory ||T^ , we consider 
the function 

N-1 

Km{X, X;t)=J2 Qm{X, t) hm{\ t) (16) 

m=0 

with the two important properties 

' dX KN{X,X;t) = N , (17) 

/ dfi KNiX,fi;t) KNil^,X;t) = KNiX,X;t) . (18) 
Jo 
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The application of a general theorem of the theory of random matrices 0,0] immediately 
yields the m-point correlation functions 

-Rm(Ai, ...,Xm]t) = det [K^lXi, Xj] t) (19) 



This expression together with Eqs. (p|, p!4| , [r6D forms the key result of this letter. For m = 1 
{m = 2) it represents the density (the two-point function) 

Ri{X;t)=Kr,{X,X;t) , (20) 
i?2(A,A;t) =i?i(A;t)i?i(A;t) (21) 
-Kjv(A,A;t)K^(A,A;t) . 

In the insulating regime t ^ 1 the fc-integration in (|l^) can be evaluated by a saddle-point 
approximation |jT9| which yields for the density of the variable x = arsinh(-\/A) a sum of 



Gaussian distributions with centers 2t(l + 2m) (m = 0, 1, . . . , — 1) and variances 2t. This 
result reflects the self averaging behavior of x in the localized regime and is well understood 
by a direct analysis of the DMPK-equation [ITsHT^ in terms of Lyapunov exponents. The 



results of a numerical evaluation of (O) for = 5 and three typical length scales {L/^ = 
0.2, 1, 10) which correspond to the metallic, crossover, and insulating regime are given in 
Fig. |I| which shows the density p{x) = sinh(2a;) i?i(sinh x, L/{2^)) versus x. Both axis have 
been rescaled with the factor L/^ for convenience. In all three cases, the effect of the level 
repulsion is easily recognized because it leads to distinct maxima and minima in the level 
density. This effect is most prominent in the localized regime (curve (a), L = 10^) where 
each level has his own Gaussian peak well separated from the other levels. In the crossover 
regime (curve (b), L = ^) the levels begin to mix but the amplitudes of the maxima are still 
rather large in comparison to the minima. Even in the so-called metallic regime (curve (c), 
L = 0.2^) the oscillations are still visible and the density is not entirely constant. On the 
other hand, a constant density is expected for this regime [jl6[. This discrepancy is a finite 



effect (note that = 5 in Fig. |I|). The usual metallic regime, / -C L <^ ^ = 2A^/, exists 
in the large A^ limit only. 

Now, we focus attention on the average conductance and its second moment. The con- 
ductance is expressed in terms of the density (|20|) as 



{g)= dX-—-J2QraiX,t)hraiX,t) . (22) 

^0 1 + A 

Since Qm{X, t) is a polynomial of degree m in the variable 1 + A, we can use the decomposition 

-1^ 



Yqrx^-^^' ^) = T+T ^ (23) 

where rm-i{X) is a polynomial of degree m — 1 that does not contribute in the integral (^2]) 
due to the biorthogonality between the Q„ and the h^- The A-integration can be done 
resulting in 



{g)=2j2 dk e-(i^+2mf+k-)t ^ ^^^^ ^nk 



2^72 + 1 
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where we have introduced the coefficient 



r(iv + i + if) r(jv + i- 



Using the expressions (120|) and (^) for the density and the two-point function, it is also 



possible to get by a straightforward calculation the second moment of the conductance 

1 ^-1 roD 

(/) = - E / c^A;e-«i+2-)'+'=')* /ttanh (f ) 

2 rn=0 •'^ 

x{2m + l) a{N,m,k) . (26) 
We note that this result can also be derived in a more direct way by the general identity 

d{g) 



dt 



-4 {{92) + f ((/) - {92))^ (27) 



that follows [0 directly from the original Fokker-Planck equation for arbitrary (3 (with 
g2 = Z]i(l + -^i)"^)- III the unitary case [3 = 2 the average ((72) drops out and {g"^) is entirely 
determined by (^ and (j^^). 

Eqs. (p^ and (p6|) are exact for all values of N and t = L/ (2^). In the insulating regime, 
they can be simplified to give 

{9) - 4 (/) ^ ^ a(iV, 0, 0) (A) -^/^ e--/^« , (28) 

where the coefficient a(iV, 0, 0) takes values between 7r/4 for iV = 1 and 1 for N 00. Eq. 
(p8|) differs considerably from the behavior of the typical conductance: exp((ln5f)) ~ e^^^/^. 

We can also consider the limit N ^ 00 for all regimes. Then, the coefficient a{N, m, k) 
tends to unity and the range of the m-sum is extended to infinity. In this case, the resulting 
expressions for the average conductance and the second moment become identical with those 
obtained from the one-dimensional non-linear cx-model (for the unitary case) by Zirnbauer 
et al. PJl0[]. (There is still a factor 2 that accounts for a different consideration of the spin 
degeneracy in the Landauer formula.) 

This exact agreement for the whole range of length scales from the metallic to insulating 
regime strongly suggests that the Fokker-Planck approach [^,0 which leads to the DMPK- 
equation is completely equivalent (in the large A^-limit) to the quasi one-dimensional micro- 
scopic models from which the one-dimensional a-model is derived P,|lT],|l2l . We have proven 
this equivalence for the ffist two conductance moments. 

This result sheds new light on the symplectic case at /3 = 4 where a similar agreement 
seems hard to believe. At the moment, exact expressions of (g) and (g^) for the DMPK- 
equation and (3 = 4 are not known but it is nevertheless clear that the typical conductance 
decreases exponentially with the length of the wire. In fact, the Lyapunov exponents are 



known for every value of /3 = 1,2,4 |13-15] and one expects for the density of the x- 



variable in the insulating regime a similar behavior as in Fig. [I| (a), with Gaussian maxima 
at the positions 2t(l + /5m). The typical conductance corresponds to the ffist maximum 
whereas the average conductance (g) = f dx p{x) cosh~^(a;) is determined by the density 
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at small x where the deviations from the Gaussian approximation are rather strong. A 
precise quantitative estimation of (g) is thus not possible within this simple picture. In Ref. 
23| it was found by an involved and approximate calculation that the average conductance 



decreases exponentially with the length of the wire for all values of /?. 

On the other hand, a simple and precise argument can be given by the exact expression 
(^). Let us assume that the average conductance approaches exponentially its "finite non- 
vanishing limit" for t — ^ cx), as suggested by the cx-model result for j3 = 4. Then d{g)/dt 
vanishes and we get from Eq. (p7|) and the inequality < g2 < immediately {g2) = {g"^) = 
thus contradicting the assumption. 

Apparently, the Fokker-Planck approach and the supersymmetric description for quasi 
one-dimensional conductors seem to be equivalent for /? = 2 and disagree in the insulating 
regime for (3 = 4. It must be emphasized that the unusual results of the cr-model approach 
PJTOd are due to a so-called "zero mode contribution" which is absolutely necessary 



to reproduce the correct behavior in the metallic regime. At the moment, it is not clear 
whether there is a particular problem with the cr-model, probably related to the large N 
limit, or whether in the Fokker-Planck approach the assumption of an isotropic diffusion 
in transfer matrix space needs revision. This assumption could in principle be wrong just 
for /3 = A. A satisfactory answer to these questions remains an outstanding and important 
problem. 

The author thanks J.-L. Pichard, C. W. J. Beenakker, and M. R. Zirnbauer for helpful 
discussions and A. Miiller-Groeling for reading the manuscript. The D.F.G. is acknowledged 
for a post-doctoral fellowship. 
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FIGURES 



FIG. 1. The density p{x) L/^oi the variable x = arsinh(\/A) versus x^/L for the cases L = 10^ 
(a), L = e (b), and L = 0.2^ (c). 

Fig. 1 
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